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ABSTRACT 

Hierarchical turbulent structure constituting a jet is considered to reproduce energy-dependent vari- 
ability in blazars, particularly, the correlation between X- and gamma-ray light curves measured in 
the TeV blazar Markarian 421. The scale- invariant filaments are featured by the ordered magnetic 
fields that involve hydromagnetic fluctuations serving as electron scatterers for diffusive shock accel- 
eration, and the spatial size scales are identified with the local maximum electron energies, which are 
reflected in the synchrotron spectral energy distribution (SED) above the near-infrared/optical break. 
The structural transition of filaments is found to be responsible for the observed change of spectral 
hysteresis. 

Subject headings: BL Lacertae objects: individual (Mrk 421) — galaxies: jets — magnetic fields — 
radiation mechanisms: nonthermal — turbulence 
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1. INTRODUCTION 

A noticeable feature associated with blazars is that the 
updated shortest var iability timescale reaches a few min - 
utes (e.g., Mrk 421: ICuil[200l iBlazejowski et al.ll2005f L 
not likely to be reconciled with the light-crossing time at 
the black hole horizon. One possible explanation for this 
fact is that small-scale structure does e xist in the parsec- 
scale jet anchored in the galactic core (jHonda fe Hondal 
2004). Indeed, in the plausible circumstance that the 
successive impingement of plasma blobs (ejected from 
the core) into the jet bulk engenders coUisionless shocks, 
electromagnetic current filamentatio n (characterized by 
the s kin depth) could be prominent (jMedvedev fc Loebl 
1999). It is known that the merging of smaller filaments 
leads eventu ally to accumulation of magnetic ener gy in 
larger scales (|Honda et alj|2000at IsiTva et al.ll2003fh 

Reflecting the self-similar (power law) characteristic 
in the inertially cascading range, the local magnetic in- 
tensity of the self-organized filaments will obey |B| ~ 
.B m (A/d)C -1 " 2 , where A and d reflect the transverse 
size scale of a filament and the maximum, respectively, 
B m = \B\\=di an d (> 1) corresponds to the filamen- 
tary turbulent spectral index. The value of d is limited 
by the transverse size of jet (or blob size; D). Then, 
it is reasonable to consider that in fluid timescales, the 
well-developed coherent fields are sure to actually meet 
hydromagnetic disturbance independent of the filamen- 
tation; that is, the turbulent hierarchy is established (see 
Fig. 1). The spectral index of the superposed fluctuations 
[denoted as f3' (> 1)] could be different from (3, and the 
correlation length scale is presumably limited by ~ A. 

At this site, the electrons bound to the local mean fields 
suffer scattering by the fluctuations, to be diffusively ac- 
celera ted by the coUisionless shocks (see lHonda fc Hondal 
l2007f) . When the acceleration and cooling efficiency de- 
pend on the spatial size scales, the local maximum ener- 
gies of accelerated electrons will be identified by A (§2), 
to be reflected in the synchrotron SED extending to the 
X-ray region. More interestingly, the spatially inhomoge- 
neous property of particle energetics is expected to cause 
the energy-dependent variability of broadband SEDs. 



Here the naive question arises whether or not this idea 
is responsible for the observed elu sive patterns of energy 
correlation of light cu rves (e.g.. | Takahashi et al.l 119961 : 
iFossati et ail l2000allbl : iBlazeiowski et all 120051 ): this is 
the original motivation of the current work. 

In the present simplistic model, light travel time effects 
would still prevent the detection of variability signatures 
on timescales shorter than D/(cS z ), where S z = S/(l + z), 
and 5, z, and c are the beaming factor of the jet, red- 
shift, and speed of light, respectively. However, if a fil- 
amented piece is isolated, having loose causal relation 
with the dynamics of a bulk region serving as a domi- 
nant emitter, an intrinsic rapid variability involved in the 
subsystem would be viable. Namely, it is inferred that 
the shorter timescale is at least potentially realized, and 
observable, unless energetic emissions from such a com- 
pact domain are crucially degraded b y synchrotron self - 
absorption and/or 77 absorption (e.g.. lAharonia n 2004). 
As is, the basic notion of the present model seems to pro- 
vide a vital clue to settle the debate as to the causality 
problem incidental to observed rapid variabilities. 

In this Letter, I demonstrate that the hierarchical 
system incorporated with the synchrotron self-Compton 
(SSC) mechanism accurately generates the time lag of 
gamma-ray flaring activity behind the X-ray, confirmed 
in the high-frequency-p eaked BLLac object Mrk 421 
(|Blazeiowski et al1l2005|) . We address that in general, 
both lag and lead can appear in X-ray interband correla- 
tions, accompanying the structural transition. The ma- 
jor transition history is argued in light of the observed 
spectral hysteresis patterns. We also work out (B m ,d), 
to provide the constraint on the field strength and D that 
should be compared with those of previous models. 

2. AN IMPROVED EMITTER MODEL WITH 
HIERARCHICAL STRUCTURE 

We consider a circumstance in which relativistic shocks 
propagate through a relativistic jet with the Lorentz fac- 
tor r, such that the shock viewed upstream (jet frame) is 
weakly to mildly relativistic. Note the relation of S ~ T. 
The overall geometry and relative size scales of the afore- 
mentioned hierarchy are sketched in Figure 1. Pro- 
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Fig. 1. — Schematic of the beamed jet including the emitting 
blobs and top view of the transverse cut of a blob region (with 
the diameter D). A number of circle-like "bubbles" symbolically 
represent the transverse section of scale-invariant filaments (with 
size A, whose maximum is d). Note that d is limited by D (for 
the values, see §4). The magnified view of a small sample domain 
illustrates the fluctuating magnetic field that scatters gyrating elec- 
trons bound to the local mean field. 

vided that the gyrating electrons trapped in the filament 
(with the size A) are resonantly scattered by the mag- 
netic fluctuations, the mean acceleration time upstream 
is approximately given by r acc ~ (3ryr g /c) [r/(r - 1)], 
where r/ = (3/2&)(A/2r g ) /3 _1 , b is the energy den- 
sity ratio of fluctuating/local mean magnetic fields (as- 
sumed to be b <C 1), r g (7, |B|) is the electron gyrora- 
dius (7 being the Lorentz factor), and r is the shock 
compression ratio. In the regime in which flares satu- 
rate, T acc will be comparable to synchrotron cooling time 
r syn(7, 1 B 1 2 ) . Balancing these timescales gives the (lo- 
cal) maximum 7 of an accelerated electron, described as 

7 *(A) = {g-^-^ gi (X/d)-^ +1 ^^ 2 y/^'\ where 
go = eB m d/(2m e c 2 ), 9l = 8Tr 2 £m 2 e c 4 /(e 3 B m ), £ = 
b(r — l)/r, and the other notations are standard. At 
7 = 7*, the electron energy distribution of the power- 
law form 71(7)^7 = K r y~ p d r y is truncated. 

For simplicity, k and £ are assumed to be spatially 
constant at the moment. Then, for (3 1 < 3 (see § 3.1), 
7* decreases as A increases (reflecting likely prolonged 
r acc and shortened T syn ), to take a minimum value at 
A = d, where the synchrotron flux density makes, up 
to the frequency of (3/4ir)(5 z j*\\ =d eB m /m e c)(= Ub), a 
dominant contribution to the F v spectrum (owing to the 
maximum magnetic intensity at the outer scale). As A 
decreases, the flux density tends to decrease, extending 
the spectral tail (due to the 7* increase). Apparently, 
this property has the F v spectrum steepening above Vh , 
whereas below the spectrum retains F v oc v~^ p ~ 1 ^ 2 . 
The frequency characterizing the spectral break can be 
expressed as v\> = 7.5 x 10 14 S 50 B^l^t; 3 / 2 d^Q Hz (for (3' = 



5/3; see § 3.1), where S50 = 8 z /50, B 



m,10 



i!6 



d/10 lb cm, £_ 4 = £ ACT 4 , and £ = b(r 



Bm/lOG, 

l)/r. 



The increase of 7* in smaller A is limited at a 
critical A c , below which escape loss dominates the 
radiative loss: the equation for the spatial limit, 
r 9 ( 7 *) ~ A/2, yields X c /d ~ ( 5r 7 2 5i) 2/(3/3+1) - By 
using this expression, one can evaluate the achiev- 
able maximum 7* value as 7*|a=a c = go(^c/d)^ +1 ^ 2 , 



for which the corresponding synchrotron cutoff fre- 
quency is v c — 8 z VQg 2 (\ c i / d)( 3 ' 3 + 1 )/ 2 = S z Pogi, where 
vq = (3/47r)(ei3 m /m e c). Also, by combining A c oc 

with F Vc oc |B(A c )|(p +1 )/V~ (p ~ 1)/2 , we read 
\n(vF v ) c /\nv c <x [(5- p)(3-(p- l)]/(3/?+l) at v ~ v c . 
More speculatively, this scaling might be reflected in 
\n(yF v )p/ lni/p for measured synchrotron flux peaks. 

3. PROPERTIES OF ENERGY-DEPENDENT SPECTRAL 
VARIABILITY AND HYSTERESIS 

3.1. X-Ray Interband Correlation 

In this context, we derive the ^-dependence of the flar- 
ing activity timescale (denoted as r). In < v < 
v Cl which typically covers the X-ray band, we have 
7* = [(47r/3)(^/<5 z )(m e c/e|B|)] 1 / 2 , which is written as 
7* (A, 1/) = (i//<5 z ^o) 1/2 (A/d)-^- 1 )/ 4 . Utilizing this, the 
expression of r syn (7*, |B| 2 ) is recast into r syn (A, ^)(= 
race ~ S.t) = (r /<5 z )(<5^ /^) 1/2 (A/d)- 3 ^- 1 )/ 4 , where 
r = 367r 2 mgC 5 /(e 4 -B„). The relation of A to v can be 
derived from the equality of 7* (A, v) = 7* (A), such that 

\( v )/d = {g^ 1 g\ /(fi '- 1 \8 z y Q /y)^-P'yW^YI^+ 1 \ 
Substituting this into r(A, v), we arrive at the result 

t{v) ~ (r /^)(^/ ( 5 z ^o) (ff - 1)/2 (ffo 3 '" 1 5r 1 ) ,T/(3 ^' ) , where 
a(p,P') = 3(/3-l)(3 -pyftafl + l)(,tf'-l)]; that is, 

rodz-W^H. (i) 

Note that a = (for j3 = 1) leads to r oc v~ x / 2 , formally 
recovering the scaling for a homogeneous model. Signif- 
icantly, the states of er < 1 and > 1 imply the appear- 
ance of the modes for which the X-ray a ctivity in a lower 
v lags that in a higher v ("soft lag"; rTakahash i et al.l 



Lgner v 

1996; Rcbillot ct al. 2006) and vice versa ('hard lag"; 
Fossati et al.l2000al), respectively, and a = 1 is the tight- 



correlation mode (Sem bav et al.l [20 021 ) . The mode flip- 
ping comes about through competing A-dependence of 
cooling and acceleration efficiency. In particular, the soft 
lag appears if 



(3 < (3 C = (4 -/3')/[3(2 -13')}. 



(2) 



The critical function (3 c {f3') is, for the key range of 
(1 <)/?' < 2, plotted in Figure 2. Note that for the 
special (3' = 2 case, a[= 3((3 - l)/(3/3 + 1)] < 1 is al- 
ways satisfied, and f3' > 2 ensures (3 > (3 C (because of 
(3 C < 1): (3' = 2 and > 2 lead to soft and hard lag, 
respectively, irrespective of the (3 value. While the in- 
dex (3 is expected to be variable (reflecting the long- 
term structural evolution of filaments; see §3.3 for de- 
tails), /?' would be a constant since a mechanism of su- 
perimposed magnetic fluctuations (§1) perhaps has uni- 
versality. One can exclude f3' > 2, which yields by 
no means soft lag, which is at odds with the observa- 
tional facts, whereupon wc can take the modified upper 
bound (indicated in Fig. 2, arrows) into account. With 
these ingredients, I conjecture the preferential appear- 
ance of (3 1 = 5/3 (the Kolmogorov-type turbulence), for 
which (3 C = 7/3. As for (3, it appears that for (3 — 2 
([Montgomery fc Liulfl979l ) and 5/3 (and given [3' = 5/3J, 
the model synchrotron spectra of vF v oc i/ - 32 and j/ - 40 
(in v\, < v < v v \ vF v oc is - 69 in v < v^) provide a 
reasonable fit to the measured ones (at flares) in the 
mid state (2002-2003) and high state (2004-2005), re- 
spectively, of Mrk 421, suggesting i/|,~2x 10 14 Hz (not 
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Fig. 2. — Phase diagram of the hierarchical turbulence. The 
function /3 c (/3') (eq. [5j; thick solid curve) divides the (J3, /3')-plane 
into the two domains that allow soft-lag (j3 < f3 c ; yellow) and 
hard-lag (/3 > /3 C ) in the X-ray band. The possible phases <E>2 
(2,5/3) and $5/3 (5/3,5/3) (large and small marks, respectively) 
involve soft-lag, as consistent with the correlations measured in 
the labeled epochs. The inset shows the power-law fit to two avail- 
able t-v data on 1 998 April 21 (for detailed data reduction, see 
Fossati ct al. 2000a), giving the indices of s = 0.11 and 0.087 (cor- 
responding colors), which determine (3(f3';s) (colored curves) for 
/3' < (s + 2)/(s + 1) (colored arrows). The characteristic curves 
lie in a domain (purple; /3'-range is indicated by thin solid arrows) 
restricted by the m easured X-ray spectral indices (and p = 1.6; 
IFossati ct al. 2000b), and give /3 = 3.1 and 2.9 (colored marks) at 
fi' = 5/3. For further explanation, see the text. 

shown in figure), as compat ible with smaller variabi lity in 
the ^-range below R band (jBiazejowski et aLl i2005). Be- 
low, we refer to these possible phases (/?,/?') = (2,5/3) 
and (5/3,5/3), which satisfy equation ([2]), as "$2" and 



5/3 



respectively. Also, we compare with the de- 



tailed data of burst decay time (in 1998; IFossati et al.l 
20003). Th e guideline is given in Figure 2: use is made 
of the translation of the measured timescale oc v s to 
f3 = [4-l3'+s(/3'-l)]/{3[2-(3'-s(/3'-l)]}. This charac- 
teristic curve for s ~ 0.1 indicates /3 ~ 3 at /?' = 5/3, and 
(3 > (3 C , as consistent with the measured hard lag. For 
these (3 = 5/3, 2, and 3, we anticipate (vF v ) p e x i/"' 8 , 
is® 9 , and v p (for p ~ 1.6; e.g., iMacomb et al.l ll995T) . 
amenable to the full X-r ay data analysis of Mrk 421 flares 
(jTramacere et al.ll2007t ). 

Let us now estimate the time lag of a soft energy band 
£l(> hvi,; h is the Planck constant) behind a hard band 
en(< hv c ) by At = t(el) — t(eyi). Here it is instructive 
to note the relation of (A c <) A(en) < K € l) (< d). Using 
the expression of Vb (eliminating £ 3 / 2 d -1 ), we obtain 



At = 1.8 6 



1/2^-3/2^-3/7 



m,10 v b,\± 



-1/14 U 
s L y T]-l hi 



(3) 



for the structural phase $2, where 1^,14 = i^/10 14 Hz, 
'7-i = [1- (eL/e H ) 1/14 ]/10-\ and e L>1 = e L /lkeV. Con- 
cerning the validity, it has been checked that, e . g., fo r 
a soft-lag episode (in 1994 May; iTakahashi et al.Hl996D . 
the measured time lag plotted against el could be more 
naturally fitted by the function ([3|) of Ar(eL,£H) (given 
en — 4 — 5kcV for ASCA), rather than the function of 

1 — (ejj / £h) ^] for the homogeneous (a = 0) 



-l/2r 



model. 

3.2. X/ Gamma-Ray Cross-Band Correlation 

The interband correlation property is reflected in 
the cross-band correlation between X- and gamma- 
rays, provided the SSC mechanism 
nant gamma-ray emitter (e.; 



and 

as a domi- 
iMaraschi et~aT1 [1992; 



iDermer &T Schlickciscr 1993). Along the heuristic (time 
independent) manner, we suppose 7 ~ 7* for scatter- 
ing electrons, and examine the correlation between an 
X-ray band e x (compared to en) and gamma-ray band 
e 7 susceptible to the inverse Comptonization of low- 
energy synchrotron photons (with €l). Here we focus 
on the feasible, Thomson regime of (el/^z)7* < rn e c 2 ; 
note that using the expression of 7*[A(el)] (§3.1), this 
range can be written as e L < S^B"^ i (C-4 2 ^r6 1 )] 2 ^ 23 
(for $2)- The Lorentz factor of the electrons that ex- 
ecute the boost of e 7 /eL = (7*) 2 is denoted as 7* = 
(e 7 /(5 2 /i^ ) 1/4 [A(e L )/d] _(,3 ~ 1)/8 - Then, simply estimat- 



ing At 7 



re 



.lilt?] ~ r(e x ) (> 0, for (3 < (3 C ) 



would be adequate for the present purpose. For con- 
venience, one may eliminate £l from 7* [transform A(eL) 
into A(e 7 )], and adopt the positive soft-lag representation 
of At X7 (= — At 7X ), so that the negative sign indicates 
gamma-ray lag. Again using v^i we find for $2 



Ar x 



where 



-1.75 



51] 



5/32 R -7/8 



■7/16 -1/32 



b,14 



{1 - 0.79e 



.//7U/16 



c 7,l 

1/14 
x,25 



?7 7X! _i days, (4) 



(e7,1^50^m,lo) 1 ^ 2; 



}/10" 



e 7 ,i 



e 7 /l TeV, and e x , 2 5 = e x /25 keV. 



The simultaneous equations ([3]) and (j4|) contain the so- 
lutions (S,B m ), for given observable quantities v\, and 
(Ar, At X7 ), as well as (cl, £h; ex> £7) inherent in detec- 
tors. 

In Figure 3 (top) for Vb.u = 2, At = lhr, 
and (£l,i, cl/eh; e x ,25) = (1, 0.2; 1), compared to Mrk 
421 (z = 0.031) meas urements (|Takahashi et al.l 119961 : 
iBIazejowski et aill2005t ) , the self-consistent numerical so- 
lution S is plotted against At X7 , given e 7 that covers a 
gam ma-ray band associated w ith the Whipple observa- 
tion (ICatanese fc W cckcs 1999T). For the allowed domain 
of S > 1 (jPiner et al.lll999n . a typical TeV range of e™ 1 ~ 
1 — 2 (susceptible to the significant variation in the mid 
state) is found to a priori restrict the domain of the ob- 
servable — At X7 to 1.4 — 2. 2 days. Surprisingly, this quan- 
titatively agrees with At X7 = — 1.8 ± 0.4 days that has 
been re vealed by multiband mo nitoring in the 2002/2003 
season ()BIazeiowski et aI1l2005[ ). In order to solidify the 
argument, the solutions for the high state with $5/3 
have also been sought. The results show that the up- 
per bound of — At X7 , at which S diverges, shifts (from 
2.2 days) to 1.7 days and the Whipple coverage e 7 ,i < 10 
restricts to — At X7 > 0.7 days; these combination yields 
— At X7 ~ 0.7— 1.7days. This is certainly compatible with 
the measured At X7 = -1.2 ± 5 days (in the 2003/200 4 
season; for the significance, see Bla zeiowski et alJ feOOS). 

3.3. Hysteresis Reversal via Structural Transition 

From the view point of activity history, it is claimed 
that, involving the fluctuations with a common [3' = 
5/3, the coherent structure, at least, in the domi- 
nant emissi on region has been i n the (3 = 2 phase 
(1994 May; ITakahashi et al.l ll996l). j3 = 3 (1998 April; 
IFossati et al.l l2000al |bj), an interm ediate phase aroun d 
(3 = 7/3 (2000 May and N ovember; ISembav et all2002ft . 
(3 = 2 (2002/2003 season; iBlazeiowski et al.ll2005D. and 
f3 = 5/3 (2003/2004 season: IBIazejowski et al.ll2005| ). to 
give rise to a hard and soft X-ray lag for (3 ^ j3 c = 7/3, re- 
spectively, and no lag for (3 = f3 c , as consistent with the 
observed correlation properties in each epoch (Fig. 2). 
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Fig. 3. — Beaming factor <5 (top) and maximum field strength B m 
(bottom) vs. the time lag of X-ray behind gamma-ray Ar x7 (< 
means gamma-ray lag) for a given gamma-ray band e 7 (labeled) as 
a parameter. The horizontal axes are common in the panels. The 
correlation is taken between e x = 25 keV and e 7 > 350 GeV com- 
pared to an RXTE band and coverage of the Whipple 10 m tele- 
scope (brick-like shaded area indicates the prohibited domains), 
respectively. The hatched bands (green) indicate the allowed do- 
mains for e 7 = 1 — 2TeV and S > 1, which cover the mea- 
sured At x7 = —1.8 ± 0.4 days (yellow) and give 8 = 10 — 92 and 
B m = 10 — 22 G at At x7 = —1.8 days (deep green bars; say, <5 = 29 
and B m = 16 G [marks] for e 7 = 1.5 TeV [dot-dashed curves]). The 
inset in the bottom panel shows the B m -S relation independent of 
e~. 

At t his juncture, the confirmed reversal betwee n clock- 
wise dTakahashi et alJll996tlRebillot et aLll2006[ ) and an- 
ticlockwise (Fossati et al. 2000b) hysteresis loops in the 
flux-spectral index plane is ascribed to the phase tran- 
sition between (3 < f3 c and > /? c , respectively. Physi- 
cally, the likely ( 3 = 2 is associated with th e prominence 
of filamentation (|Montgomerv fe Liulll979l h The smaller 
j3 = 5/3 in a high state arguably reflects strong structural 
deformation, while the larger [3 — 3 can be interpreted 
as th e dual-cascade p hase of two-dimensional turbulence 
(e.g., iKrom mcs 2002 an d references therein) transverse 
to pronounced filaments ([Honda et alJl2000bh . 

4. DISCUSSION AND CONCLUDING REMARKS 

The practical formula that constrains magnetic field 
strength is readily obtained from equation ([3]), and in 
parallel, one for $5/3 can be derived as well. We find the 
outcome that for <J> 2 and $5/3, B m must satisfy 



B <5 1/3 - 



^-^(Art-.fVi) 2 / 3 C 



>13^AT- 1 e^yVi) 2/3 G, 



4 -I/" 



(5) 



respectively, where 77* j = [1 — (cl/ch^^/IO -1 and 
At is in hours. In Figure 3 (bottom), we plot the 
self-consistent solution B m (against Ar X7 ; correspond- 
ing to 5-At X7 in Fig. 3 [top]) that obeys equation (0 



for $2 (inset) with the same parameter v alues as the 
top p anel. We see that the observed 8 > 1 (jPiner et al.l 
1999) provides the constraint for which local magnetic 
intensity (|B|) never exceeds 47 G for $ 2 (51 G for $5/3). 
Whereas a mean magnetic intensity B is not well defined 
within the present framework, the obtained scaling of 
Bm^afi 1 ' 3 — 5 seems to be reconciled with the conven- 
tional B5 1 / 3 ~ 0.1 — 1G derived from fitting a variety 
of homogen eous SSC models to the measured broadband 
SEDs (e.g.jGhisellini et alj |1998t p Tavecchio et aH ll998t 
iKrawczvnski et al.ll2001[) . 



"3/2 , _ 



7.5ia , A 6saB. 



-3/2 



b,14°50-O m ,io 



In turn, the quantity of £_ 4 dig 
(valid for f3' = 5/3; §2) is self-consistently determined. 
Making use of equation ([5]) to eliminate B m , we have d = 
2.6 x 10 16 (<55o£-4) 3/2 cm for $ 2 [2.4 x 10 16 (5 5 o£-4) 3/2 cm 
for $5/3], given the common parameter values (such as 
^6,14 = 2). To estimate d, here we call for another expres- 
sion, v c = 1.0 x 10 22 c5 50 £_ 4 Hz [independent of (/3,/3'); 
§ 2] . Using this to eliminate <5£ from the d-expression, 

we obtain the simple scaling of d — 8.2 x 10 14 z/?2i cm 
for $2 (7.5 x 10 14 ;/ 3 j 2 cm for $5/3), where u Ct 2i — 
z/ c /10 21 Hz. The size d implies the allowable minimum 
of D; e.g., v Ct 2i = 0.1 — 10 (yet involving the large obser- 
vational uncertainty) provides D\q > 10 -3 to 1 (where 
Z?i6 = £>/ 10 16 cm), as reconci l ed with the previous re- 
sults (e.g..lFossati et alj|2000bt [Krawczvnski et al.ll200ll 
IBlazeiowski et al.ll2005l) . It also turns out, from the de- 
scaling, that the range of v c ,2i < 10 2 accommodates 
^<Cl, and thereby the assumption of b <C 1 (§2). 

In addition, given an energy input into the jet, parti- 
cle density n is estimated. Assuming that electron in- 
jection operates at 7i n j <C j*\\=d (< 7*), we approx- 
imately get n ~ (k/0.6)7 1 ~?' 6 (for p = 1.6), to find 

that the steady luminosity of 10 44 ergs s _1 , which ap- 
pears to retain a dominant portion around the v^, re- 
quires n > 6 x 10 4 7 i ~ J °' 6 £' 1 ~ 6 3 -B~^ lu cm -3 (when suppos- 
ing a spherical emitting volume with the diameter of D). 
Recalling B m ^ < 5, we thus read n > 10 3 £>f 6 3 cm -3 
for ordinary 7j n j ~ 0(1); note that an upper bound 
can be given by imposing the conditions of, e.g., pair- 
plasma production (T > 1 MeV) and radial confinement 
(nT < B^/8ir), such that n < 10 s cm" 3 (suggesting 

016 > i(H). 

In conclusion, the gamma-ray lags of 1 — 2 days mea- 
sured in Mrk 421 have been nicely reproduced by the 
hierarchical turbulent model of a jet. The crucial find- 
ing is that the structural transition $2 — » $5/3 results 
in downshifting the upper bound of the observable lag 
[in a TeV (e 7) i ~ 1) band] from 2.2 to 1.7 days, in ac- 
cordance with a c l oser i nspection from 2002 to 2004 by 
IBlazeiowski et ail (|2005f) . A typical 1.8 day lag (in the 
2002/2003 season) suggests 5 = 10 - 92 and B m = 10- 
22 G (Fig. 3); the latter provides an upper limit of local 
magnetic intensity. The present model as a possible alter- 
native to the previous leptonic (e.g.. ISikora et al.lll994l : 
iBednarek fe Protherodll997tlKonopelko et al.l l2003h and 
hadronic scenarios (e.g., iMiicke fc Protheroell2001l ) will 
shed light on puzzling aspects of broadband spectral vari- 
ability. 
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